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We formulate a spherical harmonically decomposed 1 + 1 scheme to self-consistently evolve the trajectory of 
a point particle and its gravitational metric perturbation to a Schwarzschild background spacetime. Following 
the work of Moncrief , we write down an action for perturbations in space-time geometry, combine that with the 
action for a point-particle moving through this space-time, and then obtain Hamiltonian equations of motion 
for metric perturbations, the particle's coordinates, as well as their canonical momenta. Hamiltonian equations 
for the metric-perturbation and their conjugate momenta, for even and odd parities, reduces to Zerilli-Moncrief 
and Regge- Wheeler master equations with source terms, which are gauge invariant, plus auxiliary equations 
that specify gauge. Hamiltonian equations for the particle, on the other hand, now include effect of metric 
perturbations — with these new terms derived from the same interaction Hamiltonian that had lead to those well- 
known source terms. In this way, space-time geometry and particle motion can be evolved in a self-consistent 
manner, in principle in any gauge. However, the point-particle nature of our source requires regularization, 
and we outline how the Detweiler- Whiting approach can be applied. In this approach, a singular field can be 
obtained analytically using Hadamard decomposition of the Green's function and the regular field, which needs 
to be evolved numerically, is the result of subtracting the singular field from the total metric perturbation. In 
principle, any gauge that has the singular-regular field decomposition is suitable for our self-consistent scheme. 
In reality, however, this freedom is only possible if our singular field has a high enough level of smoothness. For 
a singular field with minimum quality, one can adopt the Lorenz gauge condition, which we have recast into our 
formalism: for each / and m, we have 2 wave equations to evolve odd and even parity gauge invariant quantities 
and 8 first order differential equations to fix the Lorenz gauge and determine the metric components. 



I. INTRODUCTION 

In this paper, we discuss the motion of a small com- 
pact object (idealizing a black hole or neutron star) moving 
around a much more massive, Schwarzschild black hole — 
and the gravitational waves such a system would emit. In 
gravitational-wave astrophysics, this process is often referred 
to as an Extreme Mass-Ratio Inspiral (EMRI). This problem 
has attracted a lot of attention in recent years due to the possi- 
bility of directly detecting these waves using space-based [1- 
3] and even ground-based laser interferometer gravitational- 
wave detectors [4]. In EMRIs, the small object stays for a 
long time outside of the big black hole, emitting many cycles 
of gravitational waves — even in the strong-field region very 
close to the big black hole. This wave, if detected, will pro- 
vide accurate information about the mass and the spin of big 
black hole, as well as parameters of the small object's orbit; 
one may even test whether the big black hole is indeed a Kerr 
background spacetime [5-8]. On the other hand, because it is 
the many cycles that would together lead to a detectable sig- 
nal, it will be crucial (at least for the weaker sources) to get 
very accurate models for the waveforms (e.g., at the level of 
one or a few radians for the entire waveform, which may be 
up to 10 6 cycles) in order to be able to extract them from data 
[9]. 

Because the orbiting object is much smaller in mass, one 
expects the application of black-hole perturbation theory [10- 
12], successively to higher orders in mass ratio, would be 
a viable program towards solving the EMRI problem, just 
like Post-Newtonian theory has worked for the inspiral of 
comparable-mass compact objects [13] — although a direct 
application of post-Newtonian theory to EMRIs will not be 



very efficient because one expects the most interesting waves 
to be generated when the small object moves very close to 
the big black hole, with where post-Newtonian theory breaks 
down very fast. By contrast, full numerical simulation of 
the entire EMRI spacetime would be very expensive due to 
the large separation of scales and very long integration time 
that is required to providing meaningful information regard- 
ing the evolution of the orbit. Nevertheless, the most extreme 
mass ratio achieved so far in numerical relativity simulations 
is 1 : 100 [14]. 

When computing the leading-order waveform emitted by 
a small compact object moving in a black-hole background, 
one can idealize the small object as a test particle moving 
along a geodesic, and perturbations to the black-hole space- 
time is sourced by a 5-function stress-energy tensor along that 
geodesic — with all other multipoles of the object ignored. 
However, if we would like to further calculate the evolution 
of the object via coupling to the radiation field, we need to 
consider space-time geometry near the object, which formally 
diverges if we still use the point-particle model. From this 
point of view, a regularization procedure is necessary. 

Historically there are two approaches toward regulariza- 
tion. The first involves matching the external (point-particle- 
generated perturbated black-hole) spacetime to the internal 
(perturbed Schwarzschild) spacetime of the small object — in 
a region where both are valid. This was proposed and carried 
out by Mino, Sasaki and Tanaka [15] as well as by Quinn and 
Wald [16], and later by Gralla, Pound, Poisson and others [17- 
19]. This approach, for the external spacetime of the object, 
has led to the separation of the total metric perturbation field 
into two pieces: h = /idir + ^taii; ^dir is the light-cone contribu- 
tion to metric perturbation from the point particle's 5-function 
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stress energy tensor (the Hadamard direct part) and h ta i\ is the 
contribution inside the future light cone of the particle (the tail 
part). Mino et al. [15] and Quinn and Wald [16] proved that 
the regularized radiation reaction is solely contributed by h m 
which is everywhere continuous but not necessarily smoothly 
differentiable. This approach is useful when one knows the 
Green's function of the background spacetime. 

The second method, usually called Detweiler-Whiting de- 
composition, keeps the point-particle description of the prob- 
lem, but instead separates the total metric perturabtion field 
h into a regular piece Iir and a singular piece lis [20]. The 
singular piece diverges at the particles's location but does not 
have any effect on the particle's motion. It can be obtained 
by either transforming its expression in local THZ (Thorne- 
Hartle-Zhang) coordinate system [21] to the background coor- 
dinate system or applying the Hadamard decomposition of the 
Green's function [22]. The regular field satisfies the homoge- 
neous Einstein's equation and is responsible for the geodesic 
deviation of the particle's motion in background spacetime; it 
is obtained by subtracting the singular field from the full field. 

Regularization procedures above normally provides, in a 
particular gauge, a force in terms of a geodesic orbit of a par- 
ticle. In order to obtain the evolution of the particle and the 
out-going wave, one must construct an algorithm to compute 
the force, and use it to drive one's particle trajectory away 
from geodesic motion [17]. There are two major approaches 
towards the numerical implementation. One way is the mode 
sum approach, developed by Barack and Ori [23], which de- 
composes each of the 10 metric components into spherical 
harmonics, and solve 10 coupled 1+1 wave equations for each 
(l,m). Because of the modal decomposition, metric compo- 
nent for each (l,m) is finite even at the particle's location. 
The particle equation of motion is then regularized mode-by- 
mode, by subtracting a series of regularization parameters for 
each (/,m) — these regularization parameters can be calcu- 
lated either from the singular field of Detweiler and Whiting, 
or to the direct part of Mino, Sasaki and Tanaka. This mode 
sum method has already been implemented by Warburton et 
al. [24] for Schwarzschild gravitational EMRI problem. 

In the second approach, one directly applies a 3+1 decom- 
position of spacetime, and tries to obtain the regular field 
directly — by obtaining a field Iir which is approximately 
the Detweiler-Whiting Iir near the particle, but gradually be- 
comes the full field at null infinite and near horizon. As 
shown by Vega and Detweiler [25], the field Jir satisfies a 
wave equation with out-going boudnary condition at infinity 
and horizon, but with a source that can be computed from the 
Detweiler-Whiting singular field h$. Diener and Vega [26] 
have implemented this method for a scalar particle orbiting a 
Schwarzschild black hole. In principal, this effective source 
method is also suitable for 1+1 evolution scheme. In practice, 
translating the 3+1 effective source into spherical decomposed 
form and implementing it into a working code still requires 
substantial amount of work. 

In this paper, we formulate a Hamiltonian approach towards 
the EMRI problem in Schwarzschild background, with the 
aim of providing a new angle to view this problem. We will 
only write down the equations, but not attempt to solve them 



numerically. We start by generalizing Moncrief's (already 
spherical harmonic decomposed) quadratic action of perturba- 
tive Schwarzschild spacetime to include a point particle, and 
write down a joint Hamiltonian for the particle and the spher- 
ical harmonic decomposed field. The total Hamiltonian leads 
naturally to a set of canonical equations that describe the joint 
evolution of the particle and the field. Moreover, since we are 
adopting Moncrief's formalism, gauge invariant part of the 
perturbation fields are separated out from the rest of the de- 
grees of freedom — these fields, together with lapse and shift, 
drive the rest of the fields. For each (l,m), there are 6 pairs 
of canonical quantities; 2 pairs are always gauge invariant, 
and evolve independently (but driven by the particle); among 
the other 4 pairs, 3 canonical momenta and 1 canonical co- 
ordinate correspond to the momentum and Hamiltonian con- 
straints, while the other 3 canonical coordinate and 1 canoni- 
cal momentum can be fixed by gauge choices. 

Although the Hamiltonian approach provides a rather good 
way of organizing the fields, we have not found any stand- 
alone regularization technique — and must instead adopt an 
existing one. In principle, taking the 4-D Detweiler-Whiting 
singular field h$ in any gauge, as long as their continuity sur- 
vives the differential operations required for constructing our 
canonical field quantities, they can be readily used to obtain 
effective sources for Jir, the effective regular field. However, 
the currently available singular field only allows the use of 
Lorenz gauge, which means we will have to fix that gauge, 
evolving the 8 above mentioned equations. 

This paper is organized as follows, in Sec. II we briefly re- 
view Moncrief's Hamiltonian approach for gravitational per- 
turbations of the Schwarzschild metric. After that we intro- 
duce additional terms into the action that describe the point 
particle. From this new action, we rederive the odd and 
even parity metric perturbation master equations as well as 
the point mass equations of motion in Sec. III. Note that 
for both odd and even parities, wave equations agree with 
the known master equations for Regge- Wheeler and Zerilli- 
Monrief functions with a point-particle source [27-29]. On 
the other hand, the point mass equations of motion are now 
subject to the influence of background metric and both odd 
and even parity metric perturbations. They have the physical 
meaning of the geodesic motion in the perturbed background 
spacetime. In Sec. IV, we will discuss possible ways to regu- 
larize the gauge invariant quantities and specific gauge choice, 
and hence obtain fully regularized set of equations for both 
the metric perturbations and the point mass. We conclude in 
Sec. V. 



II. REVIEW OF MONCRIEF'S HAMILTONIAN 
APPROACH 

The Amowitt, Deser, and Misner (ADM) approach to gen- 
eral relativity [48] was established more than fifty years ago. 
In this approach, the Einstein-Hilbert action is written in a 
3 + 1 format similar to a constrained Hamiltonian dynamical 
system: the spatial 3-metric components are treated as canon- 
ical coordinates, while their conjugate momenta are related 
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to components of the extrinsic curvature. The lapse and shift 
functions serve as Lagrange multipliers for the Hamiltonian 
and momentum constraints. As one varies the action with 
respect to the canonical coordinates (not including lapse and 
shift functions) and their conjugate momenta, a set of evolu- 
tion equations can be obtained. As one varies the action with 
respect to the lapse and shift functions, a set of constraint 
equations are obtained — these equations are to be satisfied 
at the initial time, and will keep being satisfied if the above- 
mentioned evolution equations are followed. This approach 
rewrites Einstein's equations as an initial-value problem; it is 
closely related to the modern development of numerical rela- 
tivity [47]. 



A. First-order perturbation of a static space-time in 3+1 form 

Starting from this section, we review Moncrief's imple- 
mentation of the ADM formalism to perturbed Schwarzschild 
spacetimes [30]. In general, for a static background space- 
time, if we take spatial slices orthogonal to the time-like 
Killing vector d t , and use integration curves of d, to identify 
points with the same spatial coordinates (on the spatial slices), 
we will have a constant 3 -metric Yu (with determinant y), van- 
ishing extrinsic curvature, a lapse N that only depends on spa- 
tial coordinates, and a vanishing shift vector Nj = 0. Here and 
henceforth in the paper, we shall use i,j,k,... = 1,2,3 to la- 
bel spatial coordinates on each slice. If we have a perturbed 
3-metric fty, their canonical conjugates p lJ , lapse perturba- 
tion N', and shift perturbation N'-, then the perturbative part of 
Einstein-Hilbert action, up to quadratic order in these pertur- 
bative quantities, can be written as: [30] 



d x 



dt 



(2.1) 



Here M" and Jtf" 1 are the Hamiltonian and momentum con- 
straints, respectively, 



M" 1 =-2p\ 



A/2 



V ? /;,' h u y -R u 



(2.2a) 
(2.2b) 



and 



M~ = r l/2 [p iJ Pij - \p 2 } + \Y x ' 2 \h i} \ k tiJ\ k 



+W' 2 



-2h,h 1} \ 



_l 7 l/2 



hliij ' Af" 



(2.3) 



Here ®/?y is the Ricci tensor associated with Jij- The co- 
variant derivative "|" here is with respect to the background 3 
metric Yij- 

The action / in Eq. (2.1) leads to a Hamiltonian system with 
constraints. In particular, variation with respect to the lapse 
function N' and shift function Nj gives rise to the constraint 
equations, 



while variations with respect to /jy and % l i gives rise to the 
evolution equations: 



dhij _ 8H T dn ij _ 8H T 



dt dx'J ' dt 
Here we have defined 



H T = / cPjc 



NJf* +N'J#" +N' i J^" 



(2.5) 



(2.6) 



B. Degrees of freedom 

Let us now count the number of degrees of freedom of this 
Hamiltonian system. Nominally, we start from 6 metric per- 
turbations, plus 6 canonical momenta, governed by 6 pairs 
(i.e., 12) equations of motion. However, we have 4 constraints 
and 4 gauge degrees of freedom at all times; so in principle we 
should be able to cut down to 4 gauge independent functions, 
or 2 pairs of canonical degrees of freedom — this is what 
Moncrief worked out explicitly for Schwarzschild. 

From a 3+1 point of view, we need to show that we indeed 
only have 4 independent data to specify for free at the initial 
time slice — and the evolution of these 4 independent data 
can already describe all solutions. For this, we note that when 
specifying the 12 initial perturbation functions, we need to 
subject them to 4 constraints, so there are 8 independent re- 
maining degrees of freedom. However, we have 3-D gauge 
within the slice, as well as an extra parameter determining the 
slicing, so we have 4 additional functions that can be used to 
reveal initial data that are actually equivalent to each other — 
with 4 truly independent degrees of freedom left. The evo- 
lution these 4 functions will be supplemented by the 4 con- 
straints and the 4 lapse/shift functions to complete the 12 func- 
tions. 

In fact, we can make the above arguments a little more 
precise. Suppose, after a canonical transformation, we can 
transform the Hamiltonian and momentum constraints to be- 
come independent canonical coordinates and momenta. More 
specifically, let us label the Hamiltonian constraint the first 
canonical coordinate Js?o, denote the conjugate momentum of 
=2o as £?o; let us then label the momentum constraints as 
•^1,2,3, and label their conjugate coordinates as =2i,2,3- In 
other words, we have 



= jr"', i= 1,2,3. 



(2.7) 



We will label the rest of the canonical coordinates ^4 5 and 
momenta ^4,5. Noting that all components of the momentum 
constraints already have vanishing Poisson brackets with each 
other, 



{jr li ,J(?' j }=0, U = 1,2,3. 



(2.8) 



we only need to make sure that the Hamiltonian constraint 
have a vanishing Poisson bracket with all components of the 
momentum constraint: 



Jf" = 0, Jf' = 



(2.4) 



{.Jf',Jf" j } =0. 7 = 1,2,3. 



(2.9) 
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Odd Parity 


Even Parity 


Lapse 






Ho 


Shift 




ho (I > 1) 


H u h*(l>l) 


3-Metric 


h 


(/ > 1), h 2 (I > 2) 


H 2 , K, h\ (/ > 1) , G (I > 2) 



Table I: List of lapse, shift and 3-metric perturbations 



This rather straightforward to show if we look at the evolution 
equation for the Hamiltonian constraint: 



d 
dt 



= }Nj + (terms not involving shift) . 

(2.10) 

Now imagine we already have vanishing constraints initially, 
then in order to guarantee vanishing constraints during the 
subsequent evolution — regardless of the shift function Nj, 

we must ensure that the Poisson bracket {ji?' ,Jf" J } vanish 
numerically. However, for linear perturbation theory, Jrff" and 
J(?" J are linear in the canonical coordinates and momenta, hjj 
and p lJ , their Poisson brackets are simply numbers (or rather, 
functions of the spatial coordinate) that do not depend on these 
perturbative fields. In this way, the numerical vanishment of 
{jf, J4?' j } in Eq. (2.10) is equivalent to Eq. (2.9) — hence 
Eq. (2.7) is always possible. 

Next, let us consider the consequence of the important 
property that once <2?q and ^1,2,3 starts from 0, they must keep 
being 0. This means their time derivatives must only contain 
themselves — which means, in the Hamiltonian, quantities 
3?o and ^1.2,3 must only multiply Qq and 2?\ 23 

• [only go and ^1,2,3] &=2t.2.3 ■ [only Q and ^1.2.3] 

(2.11) 

The absence of i2 4) 5 and ^4 5 in the above terms means that 
the derivative of £^4 5 and =^45 cannot include ^Pq or i2i,2,3- 
This, plus the fact that Js?o and ^1,2,3 vanishes, means that the 
evolution of ^4 5 and i2 4 .5 must be self contained, or: 



dt 



\ 5 ~ ^ 4 ,5&^4,5 : 



dt 



£2 4 . 5 ~ ^ 4 . 5 &^ 4 ,5 • (2.12) 



In this way, these 4 are the gauge-invariant canonical vari- 
ables. Another 4 equations are driven by the above gauge- 
invariant quantities, plus lapse and shift: 



dt 



3*0 ~ ^ 4 .5&^4 5&^'0&^l,2,3&^ / (2.13) 



-J> lj2 ,3 ~ ^ 5 &^4,5&^0&^l,2,3&^,2,3-(2-14) 

The final 4 equations are simply that the constraints vanish. 

As we shall see below, in his treatment of Schwarzschild 
perturbations, Moncrief did follow the above general prescrip- 
tion — by directly using constraints as canonical coordinate 
and momenta. Note that this structure seems rather generic, 
and does not seem to be limited to Schwarzschild or even 
static spacetimes — of course, it is another issue whether one 
can separate these into different (l,m) components. 



C. Schwarzschild Perturbations 

Let us return to perturbations of Schwarzschild. In our case, 
the background metric is 



cis- = - j 1 \dt 2 + — 



r 2 (d6 2 + sm z 0d<p z 



2/1/ 



(2.15) 

and we choose to start with constant-f slices in this 
Schwarzschild coordinate system, and adopt spatial coordi- 
nates 



(x\x 2 ,x 3 ) = {r,6J). 



(2.16) 



In this way, we have N = \J— l/g Q0 = v/T— 2M/r, and non- 
zero components of fij given by 



7n 



, 722 = r z , 733 -r 2 sin 2 9. (2.17) 



\-2M/r 



We shall use i,j,k = 1,2,3 to label these spatial coordinates, 
and write components of the metric perturbation hjj as func- 
tions of spacetime coordinates (t,r, 8,<j>) — and separate the 
angular dependence by decomposing them into scalar, vecto- 
rial and tensorial spherical harmonics [27] 



dm 



K(t,r)Uj£ + G(t,r)V£S + h 2 (t,r)W$ (2.18a) 



rim 



In i 



lC = l/fH 2 (t,r)Y' m 
h rA = h* 1 (t,r)Z l ^+h 1 (t,r)xi m 



(2.18b) 
(2.18c) 



Here we have defined / = 1 — 2M/r, and shall use A,B, ... = 
1 , 2 to label angular coordinates 



(i2 1 ,i2 2 ) = (0,(/)). 



(2.19) 



The conjugate momenta pij can be similarly decomposed, 
while the decomposition of lapse and shift perturbations N',N- 
are [27] 



=fH (t,r)Y lm 
h l t ?=Hi(t,r)Y lm 



h l %=ht[t,r)ZT+h {t,r)X l l 



lm 



rim 



(2.20a) 
(2.20b) 
(2.20c) 



Here we have used odd parity vector and tensor spherical har- 
monics Xj^tWtS, as well as even parity ones, Z^,Uj^,V^; 
their definitions can be found in [27, 28]; we have also listed 
them in Appendix A. For each (l,m), (I > 2), we have a to- 
tal of 10 independent functions characterizing 10 independent 
metric components; 4 of them are lapse and shift perturba- 
tions: Ho is lapse perturbation and Hi,liQ,ho are shift pertur- 
bations. The rest 6 functions are spatial metric perturbations: 
K,G,h 2 ,H 2) h*,hi. For even parity, there are 1 lapse pertur- 
bation function Ho, 2 shift perturbation functions H\ ,/zq and 4 
spatial metric perturbation function K,G,H 2 ,h\. For odd par- 
ity, there are no lapse perturbation, 1 shift perturbation func- 
tion Iiq and 2 spatial metric perturbation functions /ii,/i2 for 
odd parity. This counting is shown in Table I. 
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D. Odd Parity (/ > 2) 

Let us first look at odd-parity perturbations, which contain 
spatial-metric perturbations h\, h 2 [Eq. (2.18)] and shift per- 
turbation ho [Eq. (2.20)]. All odd-parity infinitesimal coor- 
dinate transformations within the spatial slice can be repre- 
sented using the odd-parity vector harmonic X A , 



Q.' A = a A + £ C lm (r, t)x' A m =£l A +C A , 

lm 



which, after applying 



Shi 



Ci\j+C 



J\> 



lead to 



8h]=C r C, 8hj = — 2C . 



(2.21) 



(2.22) 



(2.23) 



Here T 2 = is constraint; (k\ , X\ ) is the gauge-invariant sector, 
which evolves independently (once setting i 2 = 0). Gauge is 
fixed by choosing ho, which correspondingly fixes the evolu- 
tion of k 2 . (Note that % 2 is constraint and should vanish.) As 
an example, the Regge-Wheeler gauge is obtained by impos- 
ing that k 2 = 0, which requires setting 



ho = 



4A(A + 1) 



(2.30) 



These odd-parity perturbation equations determine 2 out of 
the 6 spatial 3-metric components, and 3 out of the total 10 
space time 4-metric components. 

The Regge-Wheeler function frequently used in the liter- 
ature to describe odd-parity perturbations [27-29, 41^43] is 
related to k\ by 



Moncrief defined new perturbation functions 



k\ = h\ 
which transform as 



r 



Wm = fh/r 



(2.31) 



(2.24) ^/ RW j s invariant under infinitesimal gauge transformations. 



8k x = 0, Sk 2 



-2C 



(2.25) 



In other words, k\ is invariant under infinitesimal coordinate 
transformations while k2 is vulnerable to the specific choice of 
gauge. In terms of k\ <2 , and their canonical conjugates T\ <2 , the 
the odd-parity Hamiltonian [Eq. (2.6)] can now be expressed 
as 



1 



dr 



2A(A + 1) / dr ^k x 



r 2 f 

H - 

A 

r 1 



T2- 



1 



1 



2 / drhoii 



with 



A = (/-l)(/ + 2)/2 



(2.26) 



(2.27) 



Variation of the shift function ho in the Hamiltonian gives the 
odd-parity momentum constraint equation 



T 2 =0 



(2.28) 



Equations of motion for the dynamical variables take the 
form [30] 



~di 
dk 2 

dx 2 



SH T 

: "5tT 2(A + 1) 

SH T _ 4A(A + 1) 
~~8ki~~ 7 2 
SH T fr 2 



f T 2 - 



2r 2 



2A(A + 1) 



8t 2 A(A + 1) 
SH t 



2r 2 



■2ho 



8k 2 



= 



(2.29a) 
(2.29b) 

(2.29c) 

(2.29d) 



E. Even Parity (/ > 2) 

For even parity, there are 4 spatial-metric perturbations, 
K,G,H 2l h\ [Eq. (2.18)], 1 lapse perturbation H and 2 shift 
perturbations Hi,h,Q [Eq. (2.20)]. Moncrief found it conve- 
nient to recombine K,G,H 2) h\ and define a new set of vari- 
ables qi,q 2 ,q3,q4- Like k\ for odd parity perturbation, q\ 
is invariant under infinitesimal gauge transformation whereas 
<72,<?3,<?4 are gauge dependent. The conversion between 
K,G,H 2 ,h\ and q\,q 2 ,q-i,q\ can be found in [30] and Ap- 
pendix B in this paper. In terms of the new coordinates and 
their conjugate momenta, JTi ,2,3,4, the even parity Hamiltonian 
is given by: 



Hj — 



dr\ — f%4 



r{%\ - n 2j ) - 



, M 



f 
2r 2 X 

A + l 



2A 



"3 
A + l 



-2k-$ [rA.K\ +K 2 (rA) ir ] 



4(A + 1) rA 2 
2A 2 



qiq2 



2(A + 1)A 2 



l2-qi,r) +" 



.2 A 3 



Mq 2 (q 2 ~q hr ) Mq 2 2 
2(A + l)Ar 2r r 



Hoq 2 



■Hw4+hl 



2^3 (r %t\\ r 



2.32) 



Here A is defined as 



A.2(A + 3 f)=(/-l)a + 2) + 6 f 



(2.33) 
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From this Hamiltonian, it is straight forward to obtain the fol- 
lowing canonical equations of motion: 



f A A + 1 , 

qi = - 2/7T4 + -7- 7Z 3 + — =— /A^l 
Ar A 



Xq 2 4A 2 
711 -M2 - ^ 91 



(A + 1)A 2 



M^2 



2(A + l)Ar 



2/ 2 / M\ 

92 =^2^3 - - I / - — J K4 -2(f7t4), r 



A 



#3 



M(2 g2 -gl,r) #0 

2(A + l)Ar 2 



> 2 A(A + 1) r 2 A 



2hl 



Ki = - y 92,1 



94 - 



M 
K4 



2(A + 1) 



■Hi+r 



, h* 
2 / "0 



2 / M s 

-2f{lti-7h,r)-- (/ ~ I 



M 

7T4 =- -^<?2 



(2.34a) 

(2.34b) 
(2.34c) 

2 
r 

(2.34d) 
(2.34e) 
(2.34f) 

(2.34g) 
(2.34h) 



By varying the lapse perturbation Ho and the shift perturba- 
tions H\,h.Q, it is straightforward to obtain the hamiltonian 
constraint equation 



F. Monopole and dipole perturbations 

For I < 1, the evolution of 3-metric perturbations can all 
be fixed by the constraint equations plus arbitrary choices of 
lapse and shift functions. More specifically: 

For / = 0, there are only even-parity perturbations. We 
have lapse perturbation Hq and shift perturbation Hi, plus 
metric perturbations H2 and K — while lapse perturbation 
ho* and 3-metric perturbations G and h\ all vanish due to 
the non-existence of the vector and tensor harmonics Z and 
V. However, there still exists one Hamiltonian constraint and 
one momentum constraint. We can transform (H2,K) into two 
new canonical coordinates, one of them the Hamiltonian con- 
straint, the other the canonical conjugate of the momentum 
constraint — leaving no gauge-invariant perturbation fields. 

For 1 = 1, even parity perturbation, we have lapse per- 
turbation Ho, shift perturbations H\ and h^, plus three non- 
vanishing 3- metric perturbation fields, namely H2, K and h\. 
However, there exists 1 Hamiltonian constraint and two mo- 
mentum constraints, and we can transform (H2,K,h\) into the 
Hamiltonian constraint and the canonical conjugate of the two 
momentum constraints, also leaving no gauge-invariant per- 
turbation fields. 

For 1 = 1, odd parity perturbation, we shift perturbation 
ho, and one metric perturbation field, which is hi. We also 
have one momentum constraint, therefore it means a spatial 
operation on h\ will become the canonical conjugate of the 
momentum constraint, this means we have no gauge-invariant 
perturbation field as well. 



III. 3+1 APPROACH WITH POINT MASS SOURCE 



q 2 =Q (2.35) 
as well as the even-parity momentum constraint equations 

7T 3 = 7Z4 = 0. (2.36) 



In this paper, we are interested in the joint evolution met- 
ric perturbations and the motion of a point particle. In this 
section, we will augment Moncrief's formalism with a point 
panicle. 



Note that (<7i,Ti) is the even-parity gauge-invariant sector; 
(712,^3,^4) are the gauge-dependent sector, which are deter- 
mined after the lapse Ho and shifts (H\ ,h*) are fixed. 

For example, the Even-Parity Regge- Wheeler gauge is de- 
scribed by qj = q4 = 7Z2 = 0, which requires initially setting 
<73 = q\ = %2 = 0, and keeping it true by imposing ^3 = q\ = 
JT2 = through setting the appropriate h* Q [Eq. (2.34e)], H\ 
[Eq. (2.34g)] and Ho [Eq. (2.34d)]. These even-parity pertur- 
bation equations determine the other 4 of the 6 spatial 3-metric 
perturbations, and the other 7 out of the the 10 spacetime 4- 
metric perturbations. 

Gauge invariant quantity y/zM are commonly [27-29, 41- 
43] used for even parity perturbation and it is defined by 
qi/{X + 1)/A or equivalently [42] 



A + l 



K 



r J7 



2/ 
A 



r-2 



dG 



2h\ 



(2.37) 



A. 3+1 Formulation 

Since Eq. (2.1) is the action for free metric perturbations 
alone, we need to add the action for the point particle. Using 
the prescriptions in [48], we can write: 



/,„ = m J dx 

= yVxS (3) (r-Q(f)) 
+ Lf7>(5 (3) (r-Q(f)) 



~dt 



P ) —-^{ 8 ij P i P j + m 2 ) l/2 



(3.1) 



Here, regarding quantities of space-time geometry, we have 
gij = Yij + hij the total perturbed spatial metric; jV = N + N' 
the total lapse, and ,Jfi=N\ the total shift (recall that N { = 0); 
regarding the particle, Pj (y'= 1 ,2,3) are components of the 3- 
momentum, and Q(t) represents the spatial coordinates of the 
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particle, which are, more specifically, (R(t), ©(?),<!>(?)). The 
8 function is more explicitly written as 

8W(r~Q(t)) = 8(r-R(t))8(9-©(t))8($-<t>(t)) (3.2) 

From this action, we can read off the part of the Hamil- 
tonian that involves the point particle, which includes the 
Hamiltonian of the point particle alone 



H m (Q k ,P k ) = NiQ^^/yjiQ^PiPj + m 2 (3.3) 

which describes the geodesic motion of the particle, plus the 
interaction Hamiltonian that couples the particle and metric 
perturbations, 

H iD t(Q k ,P k ,N',Nj,h ij ) 

N hijfyfrPiPm 



N JfjPlPj + m 2 



-N'jPJ 



2 y/Y j PiPj + m 2 



(3.4) 



Here we have suppressed N', N, Nj, / ; and hj/s dependence 
on Q k , for simplicity — but the reader is reminded that Q k 
enters this interaction Hamiltonian through these quantities' 
dependence on (Jr. Note that H lxA is linear in the metric per- 
turbations N', Nj and hij. 

The total hamiltonian for the combined system of point par- 
ticle plus metric perturbations is 



lm 

oddT~ 



H; 



I 111 



evenT 7 



-167r// m + 167r// ir 



(3.5) 



The 16k is actually the 2k = 2{%nG) factor in the Einstein- 
Hilbert action and we are taking the Newton's constant G to be 
unity. We have now enlarged the set of canonical coordinates 
and momenta to include [Q ,Pk). The field-alone term in H tot 
describes the free propagation of metric perturbations around 
Schwarzschild, the H m term describes the geodesic motion of 
the point particle; it is //; nt that couples the fields and the par- 
ticle together: it allows the particle's motion to drive field per- 
turbations, and field perturbations to act back onto the particle 
— with action and back-action described in a self-consistent 
way. 

Due to spherical symmetry of the background spacetime, 
we can assume that the point particle is confined within the 
equatorial plane, with = n/2, Pq = — and we only need to 
deal with R(t),<P(t). In addition, because //; nt is linear in met- 
ric perturbations, we can divide it into sum of an odd-parity 
component and an even-parity component — each component 
only involving one type of metric perturbations. 



B. Odd Parity (/ > 2) 

Odd-Parity metric perturbations are described in terms of 
odd-parity vector and tensor harmonics X^",Wj^ in Eq. (2.18) 
and Eq. (2.20). For later convenience, we denote by n and E, 



n 



W^P 8 , S: 



(3.6) 



the contractions of these harmonics with angular compo- 
nents of momentum. Plugging Odd-Parity perturbations in 
Eqs. (2.18) and (2.20) into Eq. (3.4), taking Eq. (3.6) into ac- 
count, we obtain, for each (l,m), 



odd 



, „ V7 2h l P r Z + h 2 Tl 

—hoiL 

2 y/YijP'pj + m 2 

= J dr J d 2 Q.8{r-R{t))8 (2 \Q.-£l(t)) 
-h Q Z - 



r=R(t),Cl=Q,(t) 



f2(k 1 -llC2,r+lk2)P r E+k 2 n 



^YijP'PJ+m 2 



.(3.7) 



Here we have used Q. to represent (0,0), and £2 (f ) to represent 

(0(O,*W). 

For any expression inside the definition of H°^ d , for exam- 
ple Yij, it always appears along with a 8(r — R(t)) function, 
therefore being a function of (r, 9 , <p ) instead of (R, ©, <t>) does 
not seem to make a difference. However, because derivatives 
of fields are involved, we will encounter derivatives of 8 func- 
tions in further calculations, and for a generic function 'Sir) 



9(r)8'(r-R(t)) ? &(R(t))8\r-R(t)) . 



(3.8) 



This does not indicate an ambiguity in the equations of mo- 
tion that we are ultimately going to obtain, but create inter- 
mediate steps that may differ. This requires us to be careful 
with our conventions. Here we shall use the convention that 
all terms in the integrand on the right-hand side of Eq. (3.7), 
apart from the 8 functions, only depend explicitly on (r,9,<j>), 
not on (R,®,®). 

Taking the above convention for Eq. (3.7), the new con- 
straint equation is 



T 2 = - 1 -Z8(r-R(t)) 



(3.9) 



the right-hand side is singular at the location of the point par- 
ticle but zero elsewhere. Because H mt only contains metric 
perturbations, not their conjugate momenta, only the evolu- 
tion of the momenta are affected. The evolution equations for 
T] gains an additional term of: 



<?Ti 



167T- 



8k, 



add u 

Similarly for T2, the additional term is 



\6%E^8{r-R{t)) (3.10) 



d%2 



add 



I6n n 



P° 2 



8(r-R(t))- 



1 dri 


1 






r dt 


add ^ 




add.r 



(3.11) 



It is easy to check that, up to linear order (i.e., inserting 
background geodesic equations of motion for the particle) 
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this equation is consistent with the new constraint equation 
Eq (3.9). Combining the evolution equations for k\,k2,Ti,X2 
and note y/Rw = fki/r, we can derive the following master 
equation for i//rw 



d 2 



dt 2 dr* 2 
where we have defined 

r* = r + 2Mlog(r/2M-l) 

and 



Wwfot) = Soaa&t) (3.12) 



V 



odd 



3M 
r 



(3.13) 



(3.14) 



The source term is 



determined once we know ho , h \ , h 2 or ho , k\ , k 2 and their spa- 
tial derivatives at r = R(t); here ho and k 2 are related with the 
actual gauge choice and k\ is gauge invariant. If we track back 
to the wave equation for i//rw or k\, it is easy to see that k\ 
must be discontinuous at r = R(t) in order to obtain a source 
function 8'(r — R(t)). On the other hand, the equation of mo- 
tion for dP r /dt | odd contains a term proportional to k\ r . That 
means that this equation of motion is singular because it con- 
tains S(r — R(t)), This means the full metric perturbation is 
singular at the point particle's location, and directly applying 
full metric perturbation to the particle's equations of motion 
will introduce divergence. One has to apply a regularization 
scheme before one can use these equations for computation. 
This scheme must regularize gauge-invariant quantity k\ as 
well as gauge-dependent terms h 2 ,ho, since they all enter the 
particle's equation of motion. We will discuss possible regu- 
larization methods in Sec. IV. 



>odd : 



4g/ 
(A + l)r 



fU8(r-R(t)) 
P 



2P r Z8{r-R{t)) 



(3.15) 

Here the subscript "odd" means odd parity. This source term 
agrees with the ones derived in literature [11, 27, 28] as ex- 
pected. 

On the other hand, H°^ d introduces the following additional 
terms to the rate of change of the point particle's coordinates: 



dR 
dt 



d® 
dt 



= 16,^ 

dP r 



-87T 



h 2 YlP r 

W 



1- 



(Pr) 
(P°) 2 



2 1 



(3.16a) 



odd 



dPs 



167T 



P^FZ+W] 
2/r 2 sin0(P ) 3 



167T 



-pO-[hiP r X*+h 2 W++Pt\ (3.16b) 



where components of X A ,W AB can be found in Appendix A. 
Similarly, the rate of change of the point particle's momenta 
also gains the following additional terms: 



dP r 
dt 



odd 



dt 



odd 



-16,^ 
dR 



■16^ 
d<t> 



= 16?r 



' „ , 2h Ur P'Z + h 2 , r Yl 



2P° 
1 

pT) 



l6K[2h 1 P r Z+h 2 U 

J J 

(3.17a) 

2hiP r PtX*+h 2 W**(Pt) 2 



-\6n 



2P° 



(3.17b) 



Note that such a term exists for each (/,m) with I > 2. 

From the above equations of motions it is clear that the ef- 
fect of odd-parity perturbations the test particle's motion is 



C. Even Parity (/ > 2) 



Even-Parity metric perturbations are described in terms of 



the scalar harmonics Y , vector harmonics Z A m and tensor 
harmonics Ufy and Vj^ [see Eqs. (2.18) and (2.20). For later 
convenience, we define the following quantities, 

III = U%P A P B , n 2 = Vi'^P 8 , Z' = Z l A m P A (3.18) 

which are contractions of the harmonics with angular com- 
ponents of the momentum. Even-parity metric perturbation 
fields include: 

N'=- ^f 1/2 H Q Y lm , N' r = H\Y lm , N' A = h* Z A m (3. 19a) 



i Am 



H2 
' f 



\rlr n j bn 

i , n rA 



* <jlm i Am 



■ hXZ\ 



l^A i "AB 



(3.19b) 

The (l,m) component of the even-parity Hamiltonian is 

■■ = \dr [ d 2 £l8(r~R{t))8 i2) (£l~£l(t)) 

- h* Q Z' - Hi Y lm P r + HoY'"'Po j2 
f-\P r ) 2 H 2 Y hn +2h\P r Z' 



2P° 

-r 2 (ffli+Gn 2 ) 
2P» 



(3.20) 



Here again we have used £2 to represent (0,0), and £2(f) to 
represent (©(f),4>(f)), and have defined Pq = -fP°. In addi- 
tion, Hq, Hi and Ho are the lapse and shift perturbations, they 
serve as Lagrange multipliers in the Hamiltonian; K, H 2 , h\ 
and G are 3-metric perturbations, they couple with the point- 
particle dynamical variables at its location, sourcing the in- 
teraction between the field and the test mass. The relation 
between K, H 2 , h\, G and Moncrief's qi, q 2 , ^3, qt, are shown 
in Appendix B. We recall the subtlety involving 8 function 
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and its derivative mentioned below Eq. (3.7), and note that 
all terms in the integrand of Eq. (3.20), with the exception of 
the 8 function, only depend explicitly on (r,0,<p), but not on 
(/?,©,<&). 

By varying Iiq,H\,Ho, we can obtain the new constraint 
equations: 



q 2 =-l6nP Q Y lm 8(r-R(t)) 
K A =\6nY lm P r 8{r-R{t)) 



(3.21a) 
(3.21b) 



tt 3 =87Zr z Z'8(r-R{t)) + 87zY lm [r z P r S(r~R(t))], r (3.21c) 

This means q 2 , and 7I4 are all divergent at the test particle's 
location and vanish everywhere else. 

From H^ t en [Eq. (3.20)], the evolution equation for 7Z\ gains 
the additional term of 



dlt\ 



A + 1 d% 2 



add rf dt 



dn 2 
~d7 



add 

S(r-R{t)) 



add.r 



P°f 2 

and evolution equation of K 2 gains 

4n[(P r ) 2 Yi m A + 2(P°) 2 r 2 f 2 U] 



(3.22) 



d%2 

1h 



add 



P°/(1+A)A 



8(r-R(t)) 



(A + 1)A 



2fP° 



8(r-R(t)) 



(3.23) 



Similar to the odd-parity case, the evolution equations for 
7I3, % (up to linear order) are consistent with the constraint 
equations Eq. (3.21). 

Combing the evolution equation for q i,q 2 ,7t\,n 2 and the 
constraint equations, we will find that the gauge-invariant field 
y/zM satisfies a wave equation with source term coming from 
the point particle: 



d 2 



- V 1 (r) 

r even V J 



dt 2 dr* 2 
with the potential V e ' ven given by 



W7M(r,t)=S evea (r,t) (3.24) 



4/ 
r 2 A 2 



2A 2 ( A + 1 • 
18M 2 



3M 

r 
M 
r 



(3.25) 



and the source term S e give by 



2/ 2 



[A (A - l)r 2 + (4A -9)Mr+ 15M 2 ]g" 
2 



/•A 



,5 I . a — a v.- v.- (3-26) 
(A + ljA A r 

Here the Q's are master functions describing spherical- 
harmonic decompositions of point mass stress energy tensor. 
They are defined by 



T"Y lm *dQ. = -^r 8 [r -R(t)]Y ,m * [£}(?)] , 



(3.27a) 



(3.27b) 



Q r = 



Znr 1 



T rA z'r*dQ. = 



%nP r Z' 



(3.27c) 

$=%%r 2 / T AB U^*d£l=^-8[r-R(t)}, (3.27d) 



pQ 

^'• J r T AB v , m *, _ 87tr 2 U 2 8[r-R(t)} 
! ~ A(A+T)P° ' 



A(A + 1) 



(3.27e) 



The source term in Eqs. (3.26) and (3.27) agrees with the pre- 
vious derivation of Martel and Poisson [27, 28], and here we 
have adopted their notation. 

In addition to the source term in the constraint equations 
and the field evolution equations, the particle-field interaction 
Hamiltonian also generate additional terms in the particle's 
equation of motion, which causes radiation reaction. These 
terms can be obtained by varying the interaction Hamiltonian 
with respect to point mass dynamical variables, in a similar 
manner as the odd parity case, for the canonical coordinate 
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dR 
dt 



c/4> 
dt 



167T 



dP r 



-I6n 



HiY , mf l H Y'»'P> t fP r H 2 Y»» , fh\Z' 



2P° 



pO 



pO 



\67tP r 



(P r ) 2 /fH 2 Y lm + 2h\Z'P r + r 2 (KYli + GYl 2 ) 



2(P ) 3 



H Y lm P 6 



dP* 



167T 



2r 2 P°sm 2 9 2/> 2 sin 2 0(P ) 3 
h\P r Z^ + r 2 (KU^P^+GV^P^) 



pO 



(3.28a) 

f- 1 {P r ) 2 H 2 Y Im + 2h* l P r Z' + r 2 (KUi+GU 2 ) j 

(3.28b) 



and their conjugate momentum 
dH int 



dP r 
dt 



= -\6k- 



dR 



" 1 dR ' dR 2 dR 



\6n 



dP$_ 
dt 



167T- 



dH it 



d(j> 



16n 



1 

2P° 
8?r <?P° 
(P°) 2 lm 
,dZ> 



{Pr) 2 Y , m d(fH 2 (R)) | 2PrE , d(fh\(R)) | ni «9(^(/?)) 



5/? <3fl 
f-\P r ) 2 H 2 Y hn + 2 AjP r S' + r 2 (XTIi + Gn 2 ) 

dY Im 



dR 



87T 



r\n 2 H 2 



l in 



l <?y' m 



<?n, 



<?0~ 



n 2 



d(R 2 G(R)) 
dR 



(3.28c) 



(3.28d) 



We have defined 



P°(R) = J(P r ) 2 



R 2 (l-2M/R)sm 2 1 



2M/R 
(3.29) 



This set of equations, together with the even-parity wave 
equation (3.24) and the odd-parity equations (3.12), (3.16a), 
(3.17a) and (3.17b), form a complete set of self-consistent 
evolution equations for both the point particle and the metric- 
perturbation fields. 

Similar to the odd parity case, the even parity equations 
of motion also has a divergence problem. Because the wave 
equation (3.24) for i//zm contains a source term as singu- 
lar as 8'(r — R(t)), y/zM, or q\, must be discontinuous at 
the point particle's location. According to the relation be- 
tween H 2 ,G,h\ : K and qi,q 2 ,q$,q4 shown in Appendix A, 
K contains a 8(r — R(t))-type term and H 2 even contains a 
8'(r — R(t))-type term. This means terms added to the par- 
ticle's equation of motion are all singular at the particle's lo- 
cation. As a result, one has to regularize these equations of 
motion before they can be used for actual computations. 



at these orders have been solved explicitly by Detweiler and 
Poisson [50], while their back-action to the particle's canon- 
ical equations can be obtained from expressions obtained for 
/ > 2, simply removing those terms that do not exist in these 
low I's. 



IV. REGULARIZATION OF TEST PARTICLE EQUATION 
OF MOTION 



In order to obtain regular equations of motion for the point 
particle, we must carry out a regularization procedure that ap- 
propriately removes the divergences from the metric perturba- 
tion fields. While we have not been able to find a stand-alone 
regularization procedure in the 3+1 picture, currently existing 
regularization schemes can be adapted to our formalism. In 
this section, we shall outline, but not carry out, the procedure 
with which such a regularization could be done. 



A. General Discussion 



D. Monopole and dipole perturbations 

Even though there are no gauge-invariant perturbations for 
these low-/ components, metric perturbations at these orders 
do couple to the particle. The particle's perturbation to fields 



In particular, we shall discuss how the Detweiler-Whiting 
(DW) singular-regular decomposition [20, 21] approach can 
be used to regularize our canonical equations of motion. In 
the DW approach, metric perturbation field in a small but fi- 
nite region around the point particle is decomposed into the 
sum of a regular piece (superscript "R") and a singular piece 
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(superscript "S"): 



'/IV 



h R 
/iv 



-h S 

"fly • 



(4.1) 



The singular piece h s ^ v corresponds to the deformed 
Schwarzschild solution around the small test mass as seen by 
a locally free falling observer on the background spacetime — 
it is singular as we approach the location of the point particle; 
the regular piece hR satisfies the linearized vacuum Einstein's 
Equation and is everywhere regular (although it does not sat- 
isfy the out-going boundary condition at the null infinity and 
the down-going boundary condition at the future horizon). It 
is shown that /i^ v is the appropriate singularity to remove, and 
the point particle should travel along a geodesic of the per- 
turbed spacetime that differs from the background by h R v . 

DW has shown that /z£ v can be approximated analytically in 
a local normal coordinate system built around the particle — 
such as the one introduced by Thorne and Hartle [36] and de- 
veloped to higher orders by Zhang [44] (usually referred to as 
the THZ coordinate system). Another approach towards ob- 
taining h s ^ v is through the Hadamard singular Green function, 
as carried out by Hass and Poisson [22] as well as Warburton 
et al. [24]. The computation for h s ^ v is carried out as an ex- 
pansion in the proper distance away from the particle — and 
depending on the order to which this expansion is carried out, 
the corresponding h R v will only have a finite order of smooth- 
ness. 

Among components of h s ^ v , h s n is a lapse perturbation, h s tr 
and hf A are shift perturbations, while h s rn h s rA , and h s AB are 3- 
metric perturbations. One can carry out (l,m) decompositions 
of these quantities, using the appropriate harmonics, to ob- 
tain the singular pieces of our odd-parity metric-perturbation 
fields (h^h^hf) and even-parity metric -perturbation fields 
(Hq , Hf , Kq , h f , fff ,K S i G s ) . In this way, the singular metric- 
perturbation fields come with a choice of gauge (through the 
singular pieces for lapse and shift) as well as 3-metric per- 
turbations, around the worldline of the point particle. The 
(Z,m)-decomposition coefficients of the singular metric fields 
are also referred to as regularization parameters. 

It is anticipated that the mode-decomposed versions of 
these singular metric-perturbations fields should in general be 
discontinuous or singular at the radial location of the particle 
— but it is exactly these singularities that will cancel with the 
ones we obtain for the full perturbations (i.e., /i^J, 1 ), yielding 



'/IV 



k m - 
'/xv 



'/IV 



(4.2) 



which are regular. 

More specifically in the 3+1 approach, we must first ob- 
tain the full metric (including lapse, shift and 3-metric per- 
turbations), and then subtract the singular piece — resulting 
in the regular piece. A subtlety here is the choice of gauge: 
we obtain using a particular choice of lapse and shift per- 
turbations, and the arbitrariness of the choice suggests that 
the subtraction (4.2) will yield a regular result only if the full 
metric and the singular metric are computed in gauges that are 
related to each other through a smooth transformation in the 
region near the particle. 



Note that the singular field is only defined in a region 
around the point particle — because the normal coordinate 
system (e.g., the THZ coordinate system), as well as the 
Hadamard decomposition of the Green function, is only valid 
within a distance away from the particle that is comparable to 
space-time curvature. This does not prevent us from obtaining 
a regularized set of equations of motion for the point particle, 
because for that we will only need to obtain h R v around the 
location of the particle. 

However, this has lead Vega and Detweiler (VD) [25] to 
develop a slight variant of the DW regularization approach, 
which further simplifies the regularization procedure. VD first 
assumed that we can obtain an h s ^ v that has a definition every- 
where in the spacetime, although this definition is physically 
meaningful only around the particle. They then proposed the 
application of a window function W, which is very flat around 
the location of the particle, but decays rapidly towards the 
horizon and infinity. In this way, if one defines an effective 
regular field, or effective field for short, 



tr = ,M\_ wh s 

"/IV — "/IV '' "/IV 



(4.3) 



then the effective field h R v satisfies a wave equation with a 
regular source (the full source subtracted by the result ob- 
tained by inserting /i£ v = Wh s ^ v into the wave equation), as 
well as the out-going boundary condition at the future null 
infinity and the down-going boundary condition at the future 
horizon. 

In this paper, we shall discuss how the effective-source ap- 
proach can be adapted to our 3+1 Hamiltonian formalism — 
in even- and odd-parity cases. 



B. Odd parity 

Odd parity effective fields are hoR,hiR,h2R, and their 
smoothness depends on the quality of our approximations for 
hos, his an d n 2s- We shall refer the regular field to be n-th 
order smooth if it has a smooth n-th order derivative. Right 
now, singular field is available for the regular field to have 
4-th order smoothness [49]. Let us first assume the order of 
smoothness is not an issue (e.g., assuming the singular piece 
to be available up to a rather high order), and later discuss 
options when the order of smoothness is limited. 



1. An algebraic gauge 

Out of the three metric quantities, one can construct a gauge 
invariant quantity — the Regge-Wheeler (RW) function, and 
the rest two degrees of freedom are fixed by one gauge choice 
and one constraint equation. First consider the gauge invariant 
quantity, its effective regularized piece i//rwr is given by 



f 



~ hiR+l 2( d JT--r~ h2R 



(4.4) 



while its effective singular piece y/f w (>",f) is given similarly 
by h\$, his (which are singular field components multiplied 
by the window functions). 
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The effective RW function satisfies the same wave equation 
as before [See Eq. (3.12)] but with a new source 



dt 2 



V odd (r) 



WwL(r,t) = SoddRfo*) ( 4 - 5 ) 



where the new source 5^ dd is simply the effective source, given 
by 



S ddR(r,t) = Sodd(r,t) 



d 2 
dt 2 



^2-Vo dd (r) y/Rws(r,f) 
(4.6) 

Given enough smoothness on the singular piece, a smooth 
enough iJ/rwr can be obtained by solving Eq. (4.5) and im- 
posing the out-going and down-going boundary conditions 
and infinity an the horizon. This y/RWR can then be used to 
construct the rest of the gauge-dependent fields, imposing the 
gauge condition of, for example, h\R — 0. We then obtain 
regular values for all the metric perturbation fields, as well as 
their derivatives, at the location of the point particle and will 
be able to drive its motion. To carry out this computation, we 
will need the regular field to be 2-order smooth. 

To be more specific, we can always do the coordinate trans- 
formation similar to Eq. (2.21) to shift h\R to 0. The gauge 
transformation function is given by 



C' m - -C 
r 



III! 



c 



Im 



A = XA + £ C lm X l A m = x A -r 2 dr 



lin 



Todd 
L rA 
r 2 



(4.7a) 
(4.7b) 



After the gauge transformation, according to Eq. (2.21), the 
newA' 0fi ,A' 1R! ^ s are 



h'm =0, 

h' 2R =hiR + 2r 2 / dr 



UR 



Kr =hoR-2r 2 / dr-' hlh 



(4.8a) 
(4.8b) 

(4.8c) 



Suppose the original effective field Jir is a C" function on 
the test mass's worldline, the coordinate transformation must 
be C" smooth (Eq. (4.7)) and the new effective field is C 
smooth (Eq. (4.8)). Therefore h\R — is also a viable gauge 
for evolution because it can be smoothly transformed from 
Lorenz gauge. For n > 2 the spatial derivative of the metric 
components would still be continuous. By imposing the h\R = 
algebraic gauge condition, h 2R can be immediately obtained 
by solving 



dh 



2R 



2 2r _ 
-hiR = V'rwr 
r f 



(4.9) 



As T2 is fixed by the constraint equation Eq. (3.9) and Ti can 
be obtained by solving Eq. (4.18), it is then straight forward 
to obtain Iiqr through Eq. (2.29) 



ld(h 2R 

2~ 



ns) 



t 2 



dt 



2A(A + 1) 



1 5(r 2 Ti) 



2r 2 



dr 
(4.10) 



Compared to the Lorenz gauge condition (in the follow- 
ing section), computing metric perturbations in this algebraic 
gauge is relatively easier although the effective fields are 1 
order worse in smoothness. 



2. Fixing Lorenz Gauge 

Another way to ensure the smoothness of the regular field is 
to resort to the known conclusion that if we keep the full field 
in the Lorenz gauge, the existing n = 1 singular field should be 
sufficient. This has been demonstrated by Refs. [15, 21, 37]. 

The Lorenz gauge condition, 



V^/z M v =o 



(4.11) 



where h^ v is the trace reversed metric perturbation h^ v = 
fyuv - l/2gnv h apg al3 , converts into 



! )h l -f 2 r^+Xfh 2 = (4.12) 



for the (l,m) odd-parity perturbation fields. As we break this 
into singular and regular pieces, we obtain 



dhoR , nf (M 



1 )hiR-f 2 r- 



kfh 2 R=A (4.13) 



for the effective regular field components, where A is given by 

Xfh 2S (4.14) 



dh 0S „JM 
-r— 2/ 



dt 



1 h s +f r^— 
r ) or 



Combining Eqs. (4.4), (4. 13) and (2.29), we have a set of first- 
order differential equations for Jiqr and Ji 2 r 



hoR 
h 2R 



Mn M\ 2 
M 2i M 21 





hoR 


+ 


N\ 




h 2R 




N 2 



with 



Mn =M 2 2 = 

A/ 

— + 
r 



Mi2=- 



d 2f (M 



or r \ r 



M 21 =-2, 



and 



Ni ■■ 
N 2 



A 
r 



f 



2 d _2£{M 



fr 2 
"A(A + 1) 



T2 



r \ r 

1 d{r 2 Xi) 
2r 2 dr 



(4.15) 



(4.16a) 
d ' 
2dr 
(4.16b) 

(4.16c) 



(4.17a) 
(4.17b) 



Here T 2 is fixed by the constraint equation Eq. (3.9) and Ti can 
be obtained by solving Eq. (2.29) 



d(kiR + hs) 
dt 



2(A + 1) 2A(A + l)^r 



1 d(r 2 Ti) . 



(4.18) 
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At the initial time slice t = 0, we can impose the ini- 
tial gauge condition that JiQ R (t = 0) = hoR{t = 0) = and 
Eq. (4.15) determines the gauge condition evolution later on. 
Given JiQ R ,h2R, Eq. (4.4) determines the value for h\ R and 
therefore we can obtain the full set of regularized odd-parity 
field. 



C. Even parity 

For even parity, we will follow similar procedures as the 
odd-parity case. Here we are dealing with seven effective 
field quantities: K R ,G R ,H2R,Hi R ,Ho R ,h* lR ,h^ R compared to 
three field quantities in the odd-parity case. Out of these seven 
quantities, one can construct one gauge invariant perturbation 
quantity - Zerilli-Moncrief quantity and the rest six degrees 
of freedom are fixed by three gauge conditions and three con- 
straint equations. The regular piece of the Zerilli-Moncrief 
function is given by 



WZMR 



2/ 



K R 



2/ 
A 



H2R 



dK R 



r 2 dG R 



■2h 



\R 



(4.19) 



It satisfies the following wave equation 



tions are given by 



M 



= (A + l)-£ + ( — -1 

l d{H +H 2 + 2K) 
4 dt 

0=^H -2(A + l)^i- 
r r 

1 d(H Q +H 2 -2K) 

2 rJ dr 



Hi f dH x 
2 dr 



(4.23a) 



2-™\h 2 -2JK 
dHi 

'~dt~ 



(4.23b) 



0-/ 



M 



1 Ut 



X fh-Ho , 
2 G+ ^^ ]I ' 



rdh* 
2~di 



rfidh$ 
2 dr 

(4.23c) 



Combing the above equations with Appendix A, Eq. (2.34) 
as well as the constraint Eq.(3.21), one can write down the 
evolution equation for G R ,h* lR ,7t2,Hi R ,hQ R and a combination 
of effective field functions I R = Hq R + H2R + 2K R . We also 
correspondingly define Is = Hqs + H2S + 2K$. 



' Gr ' 




' Gr ' 


h* 
n \R 




h* 
n lR 


%2 


= M' 


%2 


h 




Ir 


HlR 




H m 


h* 
. OR _ 




h* 
. OS . 



•N' 



(4.24) 



dt 2 dr* 2 



■ V even (r) 



fzMR(r,t)=S eR {r,t) (4.20) 



Non-zero components of the matrix M' in Eq. (4.24) are given 
by 



with the effective source term S e R given by 



S e R(r,t) =S e (r,t)- 



-V even {r) 



WzMs(r,t) 



dt 2 dr* 2 

(4.21) 

With outgoing wave boundary condition at spatial infinity and 
black hole horizon, one can solve the wave equation and ob- 
tain the numerical value for xjfzMR or q\ R . On the other hand, 
the effective field q2R is fixed by the constraint equation 

qiR = qi-qis = -16nP Y lm 8{r - R(t)) - q 2S (4.22) 



Similar to the odd parity case, if quality of the singular field 
is high enough, we can simply set the additional lapse Hq R 
and shifts (H\ R ,hQ R ) to zero, or to use an algebraic gauge for 
the gauge-dependent fields. Also similar to the odd parity, 
one way to limit our requirement for smoothness is to apply 
Lorenz gauge condition Eq. (4.1 1), similar to what we did for 
the odd-parity case. In this case these Lorenz gauge condi- 



2/ , _ 2 
1 ■ - 72 ' M i6 - ^2 ; 



M' 2i = 2fd r -^(f-jj,M 25 = l, M' 26 = d r - r 



(4.25a) 
2 



M' 3l = -jfd 2 , M^ 2 ^(2-r«? r ), M3 



(4.25b) 
1 



M' 45 = 


A 


M 

- r 2- 


M' 51 = 


M 


■U- 

2 ' 




V 


(AM 


Ki = 


v_ 

r 


-(3- 


M' M = 


f 
2r 





2fd r , M' 46 = 



4Mf. 



'52 



34- "J' 

(4.25c) 

4(A + 1) 
r 2 n 

(4.25d) 



(4.25e) 



(3 - 5M/r)fd r - rfdf, M> 62 = 3f 2 d r , 



(4.25f) 
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Those of N' given by 



, fn 3 fAni_ 2f$ s dG s 
1 r 2 A(A + l) rX r 2 



dt ' 



5/7 



is 



(4.26a) 



(4.26b) 



5 rA 2?1 (A + 1)A 2 



V , v dG s 2 

■2-q\,r)- ^{r— h ls ) 



2r x dr r 
M(2q 2 ~q x . r ) Hps \K (P r ) 2 Y lm 
2(A + l)Ar -^ + P»WTt) S{r - R{t)) 



{ pr\2ylm 

[ -J^o-S(r-R(t)) 



A(A + 1) 

+ [\2M-r{A + A)]q 2R + 2{r-6M)qi R , r } 

1 



4(A + 1)A 



[(A - 2X)q m . r + 2rf(q 2 R,r ~ qiR,, 



(4.26c) 

djhs dig 
dr dt 

(4.26d) 



2 3M\ ~ 2/- 

a" J ^25 ~ — K s 



2(M-r)(l+X)q lR + rf(r(2+A)-2M)(q 2R -q lR , r ) 

r 3 (l+A)A 
1.5(^05+^-^) dH lS 

2 f dr dT (4 - 26e) 




(4.26f) 

We can pick the initial condition that G R = h\ R = H\ R = 
h.Q R = I R = %2 — 0. Eq. (4.24) determines their evolution later 
on. Once G R and h* lR are known, K R and H2 R can be obtained 
using Eq. (Bl). H\ R is just I R —&2R —2K R and then all ef- 
fective fields for even parity are obtained following the above 
procedure. In reality, one may let the test particle freely evolve 
for a few cycles before turning on the radiation reaction in or- 
der for the initial junk radiation to go away. Another subtlety 
here is although G Rl h\ R J Rl H\ Rl h* aR are all regular functions, 
7T2 may actually be divergent at r — R(t ). However, as long as 
the particle trajectory does not hit the grid point (which is the 



generic case and can be guaranteed by using some numerical 
algorithm), Eq. (4.24) can still be used for harmonic gauge 
evolution. 



D. Monopole and dipole perturbations 

Although there is no wave equation in the monopole and 
dipole cases, regularization does involve these orders. We 
need to carry out steps described in the above, simply ignor- 
ing the step of solving the wave equation. This is discussed 
by Detweiler and Poisson [50], but unfortunately there is no 
known unique way to remove the singular piece of contribu- 
tion from the I = 0, 1 modes. We will leave this for future 
investigation. 



V. CONCLUSIONS AND DISCUSSIONS 

In this article we have taken a 3 + 1 Hamiltonian approach 
toward the motion of a point particle around a Schwarzschild 
black hole. For the metric perturbation fields, we have sim- 
ply adopted Moncrief's perturbative Hamiltonian (quadratic 
in these fields), and his canonical transformation to a new 
set of canonical coordinates and momenta which are either 
the Hamiltonian and momentum constraints themselves, their 
canonical conjugates, or gauge-invariant (see Sec. II B). For 
the point particle, we have inserted its own Hamiltonian 
plus an interaction Hamiltonian, with the former describing 
geodesic motion and the latter describing both (i) the particle 
souring metric perturbations and (ii) the metric perturbations 
acting back onto the particle. We have obtained these equa- 
tions of motion explicitly — decomposed into (l,m) compo- 
nents (I > 2) and even and odd parities. For (i), the equations 
we obtain agree with the previous literature, obtaining wave 
equations for gauge-invariant functions that are sourced by the 
particle. In this way, we have obtained self-consistent evolu- 
tion equations for both the particle and the metric -perturbation 
fields. In principle, depending on the lapse and shift functions 
we choose, these self-consistent equations can be written for 
any gauge. The field equations we have will be in 1+1 (f and 
r) dimensions, the gauge-invariant metric-perturbation fields 
are also decoupled from the rest of the fields. 

As can be anticipated, these set of equations are singular 
due to the use of a point particle. While we have not been able 
to find a stand-alone 3+1 approach for regularization, we have 
shown that existing regularization schemes can be adopted to 
our scheme. The most straightforward approach would be 
to use the Detweiler-Whiting's singular-regular decomposi- 
tion [21], combined with the Vega-Detweiler effective-source 
approach [25]. In the case when we have a high order ap- 
proximation of the singular field, one can (i) solve the wave 
equation for the even- and odd-parity effective gauge-invariant 
fields, and (ii) fix an algebraically simple gauge for the effec- 
tive metric, obtain all effective metric components, and cal- 
culate generalized forces acting on the particle. In this way, 
we will only have evolve one wave equation for each parity 
and each / (and all m's can be taken care of simultaneously) 
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with an effective source, and use these waves, plus regularized 
gauge-fixing terms, to drive the motion of the particle in the 
self-consistent way. However, since we have not carried out 
this computation explicitly, it is not up to which level of ap- 
proximation we shall require for the singular field — although 
it may be substantially higher than what has been required be- 
fore due to the multiple spatial derivatives used in defining the 
gauge-invariant quantities. 

In case the requirement for the singular field in (ii) turns 
out to be too high, we have proposed to replace (ii) by (ii'): 
choosing the Lorenz gauge for the full field, in which case it 
was known that the currently available approximations for the 
singular field is sufficient. In this case, in addition to the two 
wave equations, we require 2 odd-parity gauge-fixing equa- 
tions, and 6 even-parity gauge fixing equations. This will be 
equivalent to decomposing Vega-Detweiler's and Warburton 
et al. [24]'s 3+1 calculations into a 1+1 form. 

For I = 0, 1 metric perturbation fields do not have gauge- 
invariant components: a canonical transformation exists to 
transform them into either the constraints or their canonical 
conjugates. In this way, we only need step (ii) or (ii') in the 
above discussion — which we have not explicitly carried out. 
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For the identity tensor Uj^: 

U l e m e =Y lm , U$ = sin 2 GY"», ff# = t$=0. (A3) 

For the even-parity, symmetric trace-free (STF) tensor har- 
monic Vjrg'- 



lm 
66 

lm 



rim 



i 



T//m 



cos 9 sin 9 + ^ I (I + 1 ) sin 2 9 
d 2 cos 9 d 



d(j) 2 ' — dQ 1 2 



ded<j> sine d<j> 



//ill 



(A4a) 



Y lm (A4b) 



(A4c) 



And finally for the odd-parity tensor harmonic WjQ: 



W, 



lm 



ee 



sine 



cos 9 d 



d9d<j) sin0 d<j) 



/lilt 



d 2 d 

sin 9 — cos 9 ^— 

dOdtj) d<j> 



wt=wt = 



1 



1 d 2 



d 2 1 d 2 d 

sin 9 ^—^r ^—r — COS 9 -r— 

d9 2 sm9d<j) 2 d9 



rim 



(A5) 



Appendix B: Conversion of Fields for Even Parity Perturbations 

The even parity perturbation quantity K,H2,h\,G and 
1\,12,13,1A are related to each other by 



Appendix A: various tensor spherical harmonics 



Here we list the components of the vector and tensor har- 
monics X%" ,Z l j" ,Wj%,Ufy,V]$ in terms of scalar spherical 
harmonics and their partial derivatives. For the odd-parity 
vector harmonics X 1 ™, we have 



lm 



1 dY 



111! 



■ X, 



lm 



sin0 d(f) ' 
For the even-parity vector harmonic Z^ m : 



sine 



dY 



lm 



39 



jl in 



lm 



dY 



ylm 



lm 



dY 



(Al) 



(A2) 



„ _ (A-l)(A-2-2A) 12-Ii,r 

Hl ~2f7pr ~' 7l + 2(ATT) 



i\ 



rA (A + 1)A 



(<72-<2v) 



2 M 2 

f[r q 3 , r -2q 4 ]. r - -j(r q 3 . r -2q 4 ) 

— T . + (rK),-K- -(1 - — ) [r q Xr -2q 4 ] 



<J = 13 
h\ = q 4 
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